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Conparlson  is  nade  between  moasuronents  of  the  light  flux  at  the 
boundary  of  a  layer  of  polydispersc  fog-oil  snoke,  due  to  radiation 
from  a  point  source  outside  the  layer,  and  predicted  values  based  on 
tho  "Six  lilux"  nothod  of  calculating  light  scattering  by  small  particles. 
Agreement  is  found  when  the  boundary  is  absorbing, but  with  a  partially 
reflecting  boundary  (albedo  0,5)  the  experiments  show  s  greater  increase 
in  flux,  due  to  reflection,  thaji  tho  thcoiy  predicts.  Again  contrary 
to  tho  theory  tho  flux  increase  becomes  groa.tor  as  tho  angle  of 
incidence  is  reduced, 

'\ 
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LIGHT  TRANSMISSION  THROUCH  A  FO&-OIL  SMOKE  LAYER  fflTH  AH 
ADSORBING  OR  PARTIALLY  REELECTING  BOUNDARY 

By 

N.V/.  Vfootten  and  7.R.  Lane 


INTRODUCTION 


The  development  of  a  method  of  prodncint^  a  stable  r,  ell-definod 
shallow  layer  of  fog-oil  smoke  near  the  floor  of  a  ohamber(l)  has  made 
it  possible  to  compare  measurements  of  the  transmission  of  light  through 
such  a  cloud  rrith  those  predicted  from  the  theoretical  analysis  of 
radiation  transport  in  such  a  system  made  by  Chu,  Clarko  and  Churchill 
(2),  and  T/ith  some  results  of  small-scale  experiments  with  dense  disper¬ 
sions  of  latex  particles  undertaken  by  Scott  and  Churchill  (3)» 

Chu,  Clarko  and  Churchill  based  their  nathernatioal  troatment  of 
transmission  of  electromagnetic  radiation  through  a  multiple  scattering 
medium  on  a  model  in  which  the  angular  distribution  of  radiation  scattered 
by  a  single  particle  is  represented  by  six  components  -  one  forvfard,  one 
backward  and  four  perpendicular  to  the  direction  of  propagatioii.  They 
have  usod  this  "six-fliax”  representation  to  obtain  analytical  solutions 
for  certain  oases  including  that  of  a  parallel  piano  dispersion  of 
particles,  bounded  by  one  absorbing  or  partially  reflecting  surface  and 
illuminated  by  a  point  source  outside  it;  for  this  configuration  they 
have  oonputod  the  flux  incident  on  the  boundar’y  surface  as  a  function 
of  the  optical  thickness  of  the  dispersion,  the  angle  of  incidence  of  the 
radiation  and  the  boundary  roflootivity.  A  brief  outline  of  the  si.x  flux 
modal  of  multiple  scattering  is  given  in  the  Appendix  to  this  paper, 

Scott  and  Churchill  (3)  later  published  measurements  of  the  flux 
received  from  a  point  source  by  a  Stt  receiver  located  in  the  partially 
reflecting  boundary.  Their  experimental  dispersion  consisted  of  very 
uniform  (0,814 dlam)  polystyrene  latex  particles  suspended  in  an 
aqueous  solution  of  glyocl.  Heptane  was  introduced  above  the  dispersion 
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to  avoid  optiopl  disoontin\iity  at  the  interface,  end  the  source,  a  neon 
lamp  (  \  =  0,2(.27/u)  located  in  the  heptane,  wac  noved  horizontally  and 
vertically  relaxive  to  the  receiver.  The  ratio  of  the  flux  striking 
the  receiver  in  the  presenoo  and  absence  of  the  dispersion  .ms  noasured 
for  a  number  of  dispersion  depths,  houjidary  3Soflectivities,  and  souroe 
positions.  The  results  v/ere  correlated  graphically  and  found  to  bo  in 
reasonable  agreement  with  the  theoretical  solution. 

The  work  dosoribod  in  the  present  paper  vtbs  undertaken  to  ascertain 
to  what  extent  the  theorotioal  computations  are  also  supported  by 
measuromonts  made  in  a  polydisperse  oil-smoke  layer. 

SMOKE  CHAMBER  EXPERIMENTS 

A  layer  of  fog-oil  smoko,  0,5  m  deep,  was  produced  on  the  floor  of 
a  chamber  (approx.  11n  x  6ra  x  3in)  by  passing  the  output  of  a  thoraal 
generator  through  boxes  containing  crushed  solid  carbon  dioxide  as 
described  in  the  appendix  to  reference  (1).  A  light  souroe  consisting  of 
a  24  volt,  150  watt  lamp  in  water-cooled  jacket  oouldbo  placed  at  any 
height  above  the  floor  and  2;f -cosine  receivers  (consisting  of  photo¬ 
voltaic  cells,  or  photomultipliors,  covered  by  ground  opal  glass 
diffusers)  v/ure  placed  facing  upvvards  at  known  distances  from  the  ground 
zero  of  the  source  to  norsure  the  flux  falling  on  a  horizontal  surface. 

The  ra.tio  of  the  flux  in  the  presence  of  smoke  to  that  in  the  absence 
of  smoko  was  fo'jnd  from  these  measurements.  The  floor  was  normally  black 
(albedo  0,05)  but  oovild  be  overlaid  vdth  white  lining  p"per  (albedo  O.5O), 

Othor  parameters  determined  during  those  experiments  vrore  (a)  the 
optical  donsity  of  tho  smoko,  neasuimd  by  the  transmission  of  a  oollimatod- 
bcam  densitometer,  (b)  the  mass  concentration  of  tho  Emoke.  detorralncd  by 
aspiration  through  filters  of  "^^illipo^e"  material,  and  (0;  the  droplet 
size  distribution,  detcmiinod  T/ith  the  aid  of  a.  oonifuge  fitted  with  glass 
slides  coated  with  a  thin  oloophobic  film.  Those  noasuroments  \7er0  made 
at  intervals  during  each  experiment  and  were  plotted  against  tine.  The 
values  of  all  tho  p.aramotors  a.t  any  given  tine  could  bo  obtained  from 
these  graphs  and  could,  therefore,  bo  related  to  each  othor. 

RESULTS  AND  CyjiCUL'.TIONS 

The  most  convenient  ncans  of  comparing  tho  results  of  our  experinonts 
v/ith  tho  published  theoretical  values  is  in  tho  form  of  a  graph  of  the 
flux  ratio  plotted  against  the  optical  thickness  of  the  dispersion  for 
various  values  of  the  c.nglo  between  the  direction  of  tho  incident  beam 
and  the  normal  to  tho  imoeiving  surface.  Our  experimental  data  wore 
reduced  to  this  form  by  means  of  the  procoduro  outlined  below.  The 
symbols  used  are  mainly  those  given  in  (j)  to  make  for  ease  of  comparison. 
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The  angle  6  between  the  direction  of  the  Incident  boom  and  the 
normal  (Fig,l)  is  deterrained  by  the  height  of  the  source  and  the  distance 
of  the  receiver  from  ground  zero.  The  flux  ratio  P  is  the  ratio  between 
the  flux  falling  on  a  horizontal  surface  in  the  presence  of  a  dispersion 
to  that  in  the  absence  of  the  dispersion.  The  optical  thickness  T  is 
the  product  of  the  dispersion  depth  t,  the  number  N  of  particles  per  unit 
volume  and  the  effective  scattering  areao^  of  the  droplets.  The  latter 
is  determined  by  the  droplet  size  and  the “scattering  coefficient  K  ,  which 
is  the  ratio  of  the  total  flux  scattered  by  a  droplet  in  all  directions 
to  the  flux  incident  on  the  geometrical  cross  section  of  the  droplet,  K 
in  turn  depends  on  the  parameter  a  =  7rd/\  (where  d  is  the  droplet  diameter 
and  X  is  the  wavelength  of  the  light  used).  The  values  of  K  were 
obtained  from  Penndorf*s  tabulations  (5)  for  a  substance  of  the  sane 
refractive  index  as  the  fog-oil,  namely  1,50.  Prom  consideration  of  the 
spectral  distribution  of  the  radiation  and  the  spectral  response  of  the 
photocells  and  photomultipliers  used  it  v/as  calculated  that  the  value  of 
X  appropriate  to  these  experiments  was  0,60yu, 

iill  the  chamber  experiments  were  started  with  as  nearly  as  possible 
the  seme  smoke  concentration.  />.  complete  analysis  of  drop-size  distribu¬ 
tion  as  dotemined  by  the  conifuge  measurements,  and  na.BS  concentration , 
from  the  nilliporc  san^jles,  was  made  for  three  of  the  eight  experiments 
and  these  were  related  to  the  e,ttenuation  coeffioient,  o',  determined  by  the 
densitometer.  The  value  of  cr,  which  is  (dven  by  the  expression 

r  =  e  (where  t  is  the  transmission  of  a  collimated  beam  over  a 
distance  x),  has  been  used  in  all  previous  chamber  experiments  to 
oharaotcrizo  the  smoke  cloud. 

The  chiuigo  in  the  cumulative  percentage  number  distribution  was 
detemined  as  the  cloud  aged  and  this  showed  that  the  moan  droplet 
diameter  inoroasod  with  time  (Pig, 2),  The  variation  of  mass  concentra¬ 
tion  and  of  O' with  time  was  also  fo\md  and  in  Pig. 3,  the  relation  betooen 
the  two  is  given.  By  taking  the  mean  value  of  the  drop  size  in  esnh 
successive  20%  range  of  the  number  versus  size  distribution  the  relative 
masses  in  each  range  could  bo  detomined.  Those  are  given  in  Table  1, 
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Table  1 

Droplet  size  and  masa  distributions  in  fof;-oil  smoke 


0 

30 

_  _ 

60 

_ 

r - - ^ 

1 

1 

Mean 
diam,  ^) 

^  mass 

^  mass 

Average 
^  mass 

i  0-20 

0.50 

4.5 

0,62 

4.^ 

0.69 

4.5 

4,4 

20  -  40 

0.63 

8.6 

0.82 

9,4 

0.88 

9.6 

9.2 

40  -  60 

0,75 

14.5 

0.97 

15.5 

1,02 

15.6 

15.2 

60  -  80 

0,88 

23.6 

1.13 

24.5 

1,22 

24.7 

24.2 

80  -  100 

1.12 

48.8 

1.40 

46.5 

1.50 

45.6 

47.0 

(it  ±e  interesting  to  note  here  that  although  the  size  distribution 
changes  as  the  cloud  eges  the  relative  mass  in  the  20?^  ranges  remains  almost 
constant). 

From  the  total  measured  mass  and  the  relative  estimated  mass,  the  number 

n  of  droplets  of  a  given  size  range  per  unit  voltme  of  the  cloud  could  be 

determined.  The  mean  size  in  each  range  was  also  used  to  determine  cr  , 

& 

The  effective  value  of  the  optical  thickness  T  ;vas  taken  asSnc^t  for 

the  five  ranges  into  which  the  drop  sizes  were  divided.  The  effective 
value  of  0  was  taken  as  the  mean  value  found  for  these  ranges.  Hie  relation¬ 
ships  beti/een  T,  o'and  a  are  shown  in  Fig. 4. 

The  close  agreement  of  the  various  relationships  in  the  three  experiments 
(well  within  the  limits  that  could  be  expected  from  the  nature  of  the  experi¬ 
ment  and  the  methods  of  sampling  the  fog-oil  smoke)  and  the  nearly  linear 
I'elationshlp  between  the  calculated  optical  thickness  and  the  measured 
attenuation  coefficient  gave  confidence  in  the  use  of  a  as  the  oharaoter- 
izing  parameter,  Tliis  v/as  taken  as  the  reference  for  the  comparison  of  the 
flux  ratios  for  the  different  boundary  conditions  and  geometrical  configura¬ 
tions  in  the  whole  series  of  experiments, 

Ihe  flux  ratio  P  measured  in  the  experiments  when  the  source  was  above 
the  smoke  layer  were  plotted  for  given  values  of  cr  against  tan  &  and  the  best 
curve  was  fitted  to  them  (Fig»5)»  From  this  curve  the  appropriate  values 
were  plotted  against  T  (Fig,6)  giving  the  final  curve  on  which  the  comparison 
with  the  theoretical  predictions  can  be  made, 

DISCUSSION  OF  RESULTS 

For  comparison  with  a  previous  series  of  smoke -chamber  experiments  (1) 

P  was  plotted  against  ground  distance  from  the  souroe.  The  agreement  was 
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veiy  good  for  all  readings  when  the  source  was  above  the  smoke  layer  hut 
for  a  source  height  of  0,5m,  whioh  is  the  approximate  height  of  the  smoke 
layer,  the  values  of  P  were  sometimes  above  those  previously  obtained 
and  sometimes  below.  Although  the  smoke  layer  is  stabilized  at  about 
O.Sm  its  height  may  well  vary  by  a  few  centimetres  so  that  the  source 
may  sometimes  be  immersed  in  the  smoke  and  sometimes  be  outside  it;  this 
would  affect  the  fluxes  received  by  distant  receivers.  Such  slight 
variations  in  the  smoke  layer  height  v/ill  have  very  little  effect  on 
the  flux  received  when  the  source  is  located  well  above  the  layer. 

When  plotted  as  a  function  of  tan  6  the  curves  of  P  for  different 
source  heights  show  little  deviation  from  a  constant  value.  This  is  in 
accordance  with  the  findings  of  Churchill  and  his  collaborators  (3), 

COMPARISON  ?^ITH  THEORY 

Only  a  fev/  of  the  data  given  in  (3)  are  suitable  for  a  direct 
oon^arison  but  in  (4.)  additional  values  of  the  flux  ratio  as  a  fxmetion 
of  tan  &  ,  T  and  boundary  reflectivity  are  given.  The  flux  ratio  E, 
quoted  there,  hov/evor,  is  defined  as  the  ratio  of  the  flux  on  a 
horizontal  surface  in  a  dispersion  to  that  failinf  normally  onto  a 
surface  facing  the  source  in  the  absence  of  the  dispersion.  P  is  there¬ 
fore  obtained  by  dividing  E  by  the  appropriate  value  of  cos  6, 

In  Fig, 6  interpolated  and  extraoolated  values  of  P,for  a  reflecting 
boundary  of  albedo  0.50 and  at  given  values  of  T  and  tan  0,  obtained  from  the 
data  (4)  are  plotted  together  with  the  results  obtained  in  our  expeilments. 

Comparison  of  the  theoretical  values  of  P  given  in  (2)  and  (3)  for 
different  values  of  a  (in  the  former  a  =  10,0  and  in  the  latter  «=  6,0) 
indioato  that  Pisa  function  of  a ,  The  value  of  a  in  the  smoke  chamber 
experiments  varies  slightly  Tdth  abut  lies  between  4,5  and  5*5 
The  theoretical  values  quoted  are  for  a  =  6,0  so  this  V7ill  account  for 
their  being  slightly  higher  than  the  experimental  values  with  an  absorbing 
boxmdaiy.  Apart  from  this  the  smoke  chamber  results  and  theoretical 
predictions  for  an  absorbing  boundary  show  good  agreement. 

For  a  partially  reflecting  boundary  (albedo  0,50),  however,  the  smoke 
chamber  results  are  somewhat  higher  than  the  theoiy  predicts.  This  is 
contrary  to  the  trend  in  the  latex  dispersion  experiments  (3)  where  Scott 
and  Churchill  conclude  that  "the  effect  of  reflectivity  on  the  transmission 
is  essentially  independent  of  the  angle  of  incidence  but  increases  signifi¬ 
cantly  with  the  optical  thickness  of  the  dispersion".  Whilst  the  smoko 
chamber  experiments  support  the  latter  part  of  this  conclusion  con^jarison 
of  the  curves  in  fig, 6  shows  that  the  inertsase  of  flux  duo  to  the  presence 
of  a  reflecting  boundary  is  greater  at  low  values  of  tan  Q  than  at  higher 
values. 
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CONCLUSIOJ^S 

It  has  bean  shown  by  chamber  experiments  that  the  attenuation  of 
light  flux  in  a  layer  of  polydisperse  fog-oil  smoke  vdth  absorbing 
boundaries  agrees  with  the  values  prcdiotod  by  the  "Six  Flux"  thooretioal 
method  of  calculating  light  scattering  in  aerosols.  It  appears,  however, 
that  the  theory  underestimates  the  increase  of  flux  due  to  a  refleoting 
boundary,  particularly  at  narrow  angles  of  incidence. 
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Appendix  to  P.T.P. (R)34 

Outline  of  Six-Flux  Model  of  Multiple  Soatterins  by  Dispersions 

For  thin  dispersions  of  uniform  paid;ioles  or  limited  paths  cf  trans¬ 
mission  such  that  re-soattering  is  negligible,  the  transmission  and  angular 
distribution  of  scattered  radiation  can  be  described  by  a  simple  exponential 
attenuation  law  and  application  of  the  Mie  theory  of  scattering  of  elootro- 
magnetio  radiation  by  dielectric  spheres.  In  dense  or  extensive  dispersions 
such  that  the  mean  free  path  for  radiation  is  comparable  ivith  or  less  than 
the  least  dimension  of  the  system,  multiple  scattering  becomes  significant 
and  the  classical  expressions  for  scattering  are  inadequate. 


Multiple  scattering  by  a  dispersion  of  non-absorbing  particles  is 
described  mathematically  by  the  transport  equation  (1,  2;  obtained  by  an 
energy  balance  over  a  differential  volume  of  unit  cross  section  and  length 

dS^  along  the  direction  Q, 

ai(r.  n)  =:  -  No-  i  (r,  n  )  +  No-  [  i(r,  n’ )  f  (fi,  n* )  dfi‘  (1 ) 

dS„  .  ®  ®  J 

W  0 

in  T/hich 


i  (r,  )  s  the  spcoific  intensity,  defined  ps  the  energy  intercepted 

per  unit  time  by  a  ^unit  area  at  r,  per  unit  sojid  angle 
from  the  direotionll  normal  to  the  area,  where  r  is  the 
Radius  vector  fiom  the  origin  of  thr  co-ordinate  system  and 
n  is  a  unit  vector  representing  dirootion.  In  general 

i(r,  Q)  is  a  function  of  five  variables,  three  for  position 
and  tv;o  for  direction. 

f  (n  °  angular  distribution  function  for  single  scattering, 

’  which  gives  tho  fpaotion  of  energy  scattered  from  a  beam 

in  tho  dirootion  fj’  into  a  unit  solid  ongl’e  in  the 
direction  For  spherical  partiolos  tho  angulfir  dls- 
tributiqn  symmetrical  about  the  incident  direction 
and  f  (n,  n’)  reduces  to  f  (d)  ns  illustrated  in  fig.A.I , 


=  Ni-mber  of  particles  per  unit  volume  of  the  dispersion. 

=  effective  cross  section  of  a  particle  for  scattering, 

Tho  scattering  cross  section  and  t?io  density  of  particles  can  be 
combined  to  give  the  mean  free  path  for  scattering  The  cross 

section,  the  moan  free  path  and  the  angular  distrlbution“f unction  are 
functions  of  the  dimensions  and  refractive  index  of  the  particles  and  of 
the  tvavolength  of  tho  radiation  (Mie  theory). 


N 


cr 

s 


Exact  solutions  of  the  transport  equation  have  been  worked  out  for 
only  a  few  oases  and  for  highly  idealized  conditions  -  isotropic  scattering 
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in  dispersions  of  infinite  thickness  -  but  useful  approximete  solutions 
can  be  obtained  if  the  angular  distribution  from  a  single  scattering  is 
represented  by  a  number  of  discrete  components.  In  the  six  flux  model 
of  Chu  and  Churchill  (3),  illustrated  in  Pig.AI ,  the  angular  distribution 
function  for  single  scntterijig  by  spherical  particles,  f(o),  is  resolved 
into  a  forward-scattering  component  f,  a  backward  scattering  component 
b,  and  four  equal  si deways -scattering  components  s,  defined  ns 

f  Vs 

f  =  27?  /  f(0)  cosV  sin  d  de 

0 

jit 

b  =  277/  f(d)  oos®0  sin  6  dd 

V 

and  s=x(l-f-b) 

(f(d)  is  the  fraction  of  radiation  scattered  into  a  unit  solid  angle 
in  the  direction  (  d),  see  Tig.Al), 

The  specific  intensity  v/itliin  a  dispersion  can  thus  be  represented 

by  six  discrete  components,  . .  Eg,  in  the  positive  and  negative 

direction  in  three  orthogonal  directions.  The  representation  for  the 
case  of  obliquely  incident  radiation  from  a  source  outside  a  plane-parallel 
dispersion  is  shov/n  in  Fig,A2,  Use  of  this  six  flux  model  permits  reduction 
of  the  transport  equation  to  a  set  of  six  simultaneous  differential 
equations,  whoso  solution  yields  the  six  components  of  the  flux  at  any 
point. 


The  solutions  obtained  by  Chu  and  Churchill  for  the  case  of  obliquely 
incident  radiation  from  a  source  outside  a  plane-parallel  dispersion 
having  a  partially-reflecting  boundary  can  bo  expressed  functionally  as 

h  =  sT,  R)  i  =  1,  2,  3,  4,  5,  6 

v/here  Z  =  No*  z  =  normal  distance  (z)  into  the  disporsJ.on  expressed  in 
moan  free  paths 

T  s  =  thickness  (t)  of  the  dispersion  expressed  in  mean 
free  paths. 

0  =  azimuth  angle  of  the  incident  radiation 

R  s  reflectivity  of  the  boundary. 

The  analytical  solution  is  extremely  complex  and  computations  are 
feasible  only  V7ith  the  aid  of  an  automatic  computer.  Chu,  Clarke  and 
Churchill  (4)  have  undertaken  such  computations  for  a  number  of  partioxilar 
cases  and  have  published  graphical  correlations  illustrating  the  effects 
of  the  various  variables  on  the  ratio  of  the  flux  at  any  point  in  the 
presence  and  absence  of  a  dispersion  of  particles  of  given  optical 
properties, 
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GEOMETRICAL  REPRESENTATION  OF  TRANSMISSION 
THROUGH  A  SCATTERING  LAYER. 


FIG.L 
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FIG. 2.  SIZE  DISTRIBUTION  OF  DROPLETS  IN  FOG -OIL  SMOKE. 
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RELATIONSHIPS  BETWEEN  T.  S’ AND 
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ESTIMATION  OF  MEAN  VALUES  OF  P  IN  TERM 
OF  Cr  AND  TAN  0. 
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FIG.AI.  ANGULAR  DISTRIBUTION  OF  RAD!.»TI0M  FROM  A  SINGLE 

scatteMgI  particle  diahete~r  comparable  with  WAVE- 

~LENGTt}:i(<x)  CONTINUOUS  DISTRIEUTION 
[b]  SiX  FLUX  REPRESENTATION 


FIG.A2.  SIX-FLUX  REPRESENTATION  OF  THE  INTENSITY  IN  A  PLANE -PARALLEL 
DISPERSION  HAVING  A  REFLECTING  BOUNDARY  WITH 
ObLlQUELY  INCIDENT  RADIATION  [  Eg  AND  Es  ARE  PERPENDICULAR 
TO  THE  PLANE  OF  THE  PAPER.’l 
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